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DTIC  QUALITY  INSPECTED  I 


Abstract 


The  dependence  of  the  second  harmonic  generation  (SHG)  signal  of  an  amorphous 
polymer  system  containing  nonlinear  optical  (NLO)  chromophores  on  the  polarizing  electric 
field  is  investigated.  It  is  shown  that  under  a  weak  poling  field  condition,  the  SHG  intensity 
is  proportional  to  Ep2,  provided  that  the  polarizing  electric  field  Ep  is  greater  than  a 
threshold  field,  Es.  Above  Es,  the  macroscopic  polar  orientational  order  is  established.  The 
lengthening  of  the  relaxation  time  for  polar  orientational  order  with  increasing  E  is  also 

Mr 

observed.  Theoretical  analysis  shows  that  the  relaxation  time  is  proportional  to  the  polar 
order  initially  induced  by  the  polarizing  field.  Experiments  have  been  carried  out  to  confirm 
the  theoretical  prediction. 


Introduction 


Amorphous  polymers  incorporated  with  chromophores  having  large  second  order 
hyperpolarizability  have  emerged  as  a  new  class  of  electro-optic  materials.  The  nonlinear 
optical  (NLO)  effect  responsible  for  the  electro-optic  response  and  second  harmonic 
generation  can  be  induced  by  an  external  electric  field,  by  which  the  macroscopic  center  of 
inversion  symmetry  is  destroyed.  The  external  electric  field  polarizes  the  sample,  and  creates 
polar  order  in  the  amorphous  polymer  system.  The  electric  field  is  usually  applied  at 
temperatures  higher  than  the  glass  transition  temperature  (Tg)  of  the  polymer,  and  removed 
at  temperature  below  Tg  when  the  induced  polar  order  is  frozen  in. 

Two  types  of  polymeric  second  order  NLO  materials  have  been  prepared:1 
X)  guest/host  systems  with  NLO  chromophores  dispersed  in  the  amoiphous  polymers; 

2)  NLO  chromophores  covalently  linked  to  the  polymer  backbone.  Due  to  chemical  bonding, 
type  2  materials  build  in  significant  impediment  to  relaxation  of  the  polar  order  induced  by 
the  electric  field.  Further  stability  in  the  orientational  polar  order  can  be  accomplished 
through  structure-enhancing  hydrogen  bonding  or  chemical  cross-linking  of  the  NLO 
chromophores  to  neighboring  polymer  chains.2  However,  due  to  ease  in  varying  the 
relaxation  rate  and  in  changing  continuously  the  NLO  chromophore  concentration  before 
the  phase  separation  sets  in,  type  1  materials  provide  as  a  useful  system  for  investigating  the 
relaxational  dynamics  of  the  second  order  NLO  effect. 

The  effectiveness  of  poling  will  determine  the  subsequent  electro-optical  activity  of 
the  sampie.  To  prepare  a  stable  NLO  polymeric  material  for  device  applications,  one  needs 
to  understand  effects  associated  with  thermal  annealing  and  electric  field  poling.  Isothermal 
annealing  has  an  effect  on  lengthening  the  relaxation  time  of  the  orientational  polar  order 
parameter  (POP)  and  thereby  stabilizing  the  second  harmonic  generation  (SHG)  response.3 


The  lengthening  of  the  relaxation  time  is  due  to  densification  of  the  material  below  the  glass 
transition  temperature,  in  connection  with  physical  ageing  of  the  amorphous  polymer  in  the 
nonequilibrium  state.  We  are  concerned  in  this  paper  with  the  effect  of  electric  field  poling. 
We  have  investigated  the  electric  field  dependence  of  the  orientational  relaxation  time  of 
the  NLO  chromophores  under  an  isothermal  condition.  In  addition,  the  SHG  intensity  is 
also  measured  as  a  function  of  the  strength  of  the  poling  electric  field.  We  show  that  both 
the  intensity  and  the  relaxation  time  increases  with  increasing  the  poling  field  strength.  The 
increase  in  the  relaxation  time  with  increasing  the  poling  field  strength  clearly  indicates  that 
the  decay  of  the  SHG  response  is  closely  associated  with  the  induced  orientational  polar 
order,  thereby  suggesting  that  the  decay  of  the  second  order  NLO  response  can  be  impeded 
by  poling  the  amorphous  NLO  polymer  with  a  high  electric  field  strength.  We  have  provided 
a  theoretical  analysis  to  substantiate  the  electric  field  dependence  effect. 

Experimental 

Appropriate  amounts  of  NLO  chromophores  NAS  (4-amino-4'-nitro-trans  stilbene) 
and  PMMA  (polymethylmethacrylate)  were  dissolved  in  chloroform  to  form  solutions  with 
different  NAS  concentrations.  After  filtering  each  solution  to  remove  undissolved 
particulates,  polymer  films  were  prepared  by  spin  coating  the  solution  on  ITO  (indium  tin 
oxide)  coated  soda  lime  glass  slides.  The  glass  slides  containing  polymer  films  were  first 
placed  in  a  vacuum  oven  at  50°C  for  over  24  hours,  and  then  baked  for  another  24  hrs.  at 
higher  temperature  (slightly  above  Tg)  to  remove  the  solvent  used  in  spin-coating.  The 
concentration  of  the  NAS  chromophore  was  monitored  by  UV-visible  absorption 
spectroscopy  and  appeared  to  suffer  little  decrease  by  the  soft  and  hard  baking  processes. 
Two  identical  glass  slides  coated  with  NLO  polymer  films  were  then  placed  in  the  vacuum 
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oven  at  120°C  to  fuse  in  a  sandwich  configuration  for  future  contact  electrode  poling. 


V 


The  second  harmonic  generation  SHG  intensity  was  measured  by  using  a  setup 
similar  to  that  reported  previously.4  It  consisted  of  a  Nd:YAG  laser  (Spectra  Physics  GCR- 
11,  A.  =  1.06  nm,  Q-switched  at  10  Hz)  used  as  an  excitation  source.  The  polymer 
film/electrode  assembly  was  placed  on  a  goniometer  stage  located  inside  a  temperature 
controlled  oven.  The  film  temperature  was  accurate  to  ±  0.2°C,  controlled  by  a  temperature 
controller.  The  incidence  angle  of  the  fundamental  beam  at  1.06  jim  was  60°;  the 
transmitted  SHG  signal  at  532  nm  was  selected  by  a  monochromator  and  detected  by  a 
photomultiplier,  followed  by  a  preamplifier  and  averaged  by  a  boxcar  integrator. 

The  glass  transition  temperatures  (Tg)  of  the  NAS/PMMA  samples  were  determined 
by  using  a  DSC  (Perkin  Elmers,  Delta  series).  Shown  in  Fig.  1  is  the  variation  of  T  with 
the  NAS  concentration.  Over  the  0-10  wt  %  concentration,  the  glass  transition  temperature 
appears  to  decrease  linearly  with  increasing  NAS  concentration.  The  refractive  index  and 
thickness  of  the  polymer  film  were  measured  by  using  a  prism  coupler  (Metricon).  The 
prism  coupler  was  operated  in  accordance  with  the  optical  waveguide  principle  where  the 
polymer  film  serves  as  the  propagation  layer  in  the  slab  waveguide  configuration. 

Results  and  Discussion 

1.  SHG  Intensity  and  Poling  Field  Strength 

The  SHG  signal  of  the  NLO  film  under  the  isothermal  condition  was  measured  as 
a  function  of  the  strength  of  the  electric  field  employed  to  polarize  (or  pole)  the  sample. 
Shown  in  Fig.  2  is  the  plot  of  the  amplitude  of  the  SHG  signal  for  the  5  wt%  NAS  in 
PMMA  sample  at  110°C  as  a  function  of  the  square  of  the  poling  field  strength  (Ep2).  One 
notes  that  there  appears  a  threshold  voltage  of  about  0.9xl06  V/cm,  below  which  the  signal 
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only  increases  slightly  with  increasing  poling  field  strength.  Above  this  voltage,  the  SHG 
signal  increases  more  rapidly  and  it  clearly  increases  linearly  with  Ep2. 

The  existence  of  the  threshold  poling  field  (Es)  has  not  heretofore  noticed,  and  it 
indicates  that  the  poling  process  proceeds  in  two  separate  stages.  The  poling  field  which 
provides  surface  charges  first  aligns  the  molecules  adjacent  to  the  electrodes,  and  the  field 
from  the  aligned  dipoles  and  additional  surface  charges  further  orient  the  remaining  NLO 
chromophores.  Because  of  thermal  motion,  the  alignment  takes  place  in  competition  with 


the  intermolecular  potentials  that  act  among  the  NLO  dipoles  and  between  the  dipoles  with 
the  polymer  host.  It  appears  that  below  Es,  the  poling  field  provides  only  enough  surface 
charges  to  align  the  molecules  near  the  surface  and  very  little  macroscopic  order  is  induced, 
hence  giving  rise  to  a  small  SHG  intensity.  Above  E&,  the  surface  charges  provided  by  the 


poling  voltage  as  well  as  the  oriented  dipoles  combine  to  overcome  the  orientation¬ 
randomizing  molecular  motion  and  induce  macroscopic  nonlinear  polarization  for  SHG. 

The  SHG  intensity  is  proportional  to  the  square  of  the  second  order  macroscopic 
susceptibility  x^2'  (~2o>;  g>,  to).1  Here  o  is  the  frequency  of  the  fundamental  laser  beam  (the 
frequency  of  the  YAG  laser  oscillating  at  1.06  jum  in  the  present  case).  Two  elements  of 
are  involved  in  the  electric  field  poled  amorphous  NLO  polymer.  They  are  xJV  and 
Xzxx2\  where  z  refers  to  the  direction  of  the  poling  electric  field,  i.e.  the  direction 
perpendicular  to  the  polymer  film  surface  on  the  ITO  substrate;  the  x  direction  is  on  the 


surface  of  the  film.  By  using  different  polarizations  of  the  exciting  laser  beam  the  values  of 
X-rrP^  and  XzxJ^  can  be  determined.4 

These  x^  elements  are  related  to  the  orientational  polar  order  parameters  L33  and 


x  J2)  -X33(2>  - 


P  ¥LX 
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and 


*J2)  -*31(2)  -  P  F  hi  (2) 

where  p  is  the  number  density  of  NLO  chromophores;  F  is  equal  to  t- 3  hv  Here  f^j  and 
f2&)  are  local  field  factors  at  frequency  o  and  2ti>,  respectively.  In  the  Lorentz  model,  the 
local  field  factors  are  given  by  fv  =  (nv2  +  2)/3,  where  r^  is  the  refractive  index  at 
frequency  v.  0  is  the  dominant  component  of  the  molecular  second  order  hyperpolarizability 
tensor.  In  writing  Eqs  (1)  and  (2),  one  neglects  the  contribution  from  other  minor 
components  of  the  hyperpolarizability  tensor. 

The  pertinent  orientational  polar  order  parameters  (POP)  L33  and  are  given  by5 

L33  (t)  =  <cos301  (t)>E  =  ~  {3<P1  (cos01(t))>E  +  2  <P3(cos01(t))>E}  (3) 

L31  (0  =  (<cos0x  (t)>E  -  <cos30(t)>E)/2  = 

|  {  <P1(cos01(t))>E-  <P3(cos01(t))>E}  (4) 

where  0X  refers  to  the  polar  angle  of  the  principal  axis  of  representative  NLO  chromophore 
1  with  respect  to  the  external  poling  field.  Here  P^cos  0X)  and  P3(cos01)  are  the  Legendre 
polymomal  of  order  1  and  3,  respectively.  The  time  dependence  in  L33  and  L^  is  caused 
by  molecular  reorientation  which  leads  to  relaxation  of  POP.  The  angular  brackets  with  a 
subscript  E  indicated  that  the  ensemble  average  is  taken  on  the  system  in  the  presence  of 
the  poling  field  E. 

As  an  approximation,  POP  parameters  L33  and  L^1  are  often  described  in  terms  of 
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Langevin  functions  of  order  3  and  l.6  Such  a  simplistic  description  is  only  valid  at  very  low 
NLO  chromophore  concentration  when  the  intermolecular  interaction  between  dipoles  can 
be  neglected.  In  reality  these  are  complex  functions  of  the  external  poling  field  and  the 
anisotropic  (angular  dependent)  intermolecular  interaction  potential.5  The  dynamics  of 
molecular  motion  renders  the  time  dependence  of  order  parameters  and  gives  rise  to  the 
relaxation  of  L33  and  L^.  However,  as  shown  in  ref.  5,  under  a  weak  poling  field  condition, 
<P3(cos01)>  is  not  excited.  In  the  next  section,  we  show  by  linear  response  theory  that  the 
time  dependence  of  POP  in  the  weak  poling  field  condition  is  given  by 

Lj^t)  =  L33(t)/3  =  -  [<cos0j(t)  cos01>  +  (N-l)  ^osO^t)  cos02>]  (5) 


where  a  =  fQ  jjEp/kT  is  the  poling  constant.  f0  is  the  Onsager  local  field  factor  associated 
with  the  poling  field,  and  m  is  the  ground  state  dipole  moment  of  the  chromophore.  0X  and 
02  correspond  to  polar  angles  for  representative  molecules  1  and  2,  respectively.  N  is  the 
total  of  number  of  NLO  chromophores  which  have  been  polarized  by  the  external  field.  The 
angular  brackets  without  the  subscript  E  in  Eq.  (5)  correspond  to  an  average  over  an 
equilibrium  distribution  function  in  the  absence  of  the  poling  field. 

For  t  =  0,  Equation  (5)  reduces  to  the  initially  induced  orientational  order 
parameter,  given  by 


L^O)  =  Lj3(0)/3  =  —  [  <cos201>  +  (N-l)  <cos0!  cos02>  ] 


a 

15 


Sk 


(6) 


Equation  (6)  was  obtained  in  ref.  5  by  using  the  method  of  equilibrium  statistical  mechanics. 
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Here  gK  is  the  Kirkwood  g-factor  given  by 


§K  =  1  + 


(N- l)<cos  0,  cos02> 
<COS20j> 


=  1  +  p  G/ 


(7) 


where  GA  is  the  cluster  integral  defined  previously. 

As  one  sees  in  Eq.  (6),  the  weak  poling  field  condition  the  POP  is  proportional  to 
the  poling  constant  a,  which  is  proportional  to  Ep;  one  thus  expects  that  the  SHG  intensity 
is  proportional  to  Ep2,  provided  that  the  poling  field  is  larger  than  the  threshold  field  Es, 
as  shown  in  Fig.  2. 

By  increasing  the  temperature,  we  have  found  that  the  SHG  intensity  first  decreases 
in  accordance  with  the  1/T  dependence  in  the  poling  constant  a.  However,  as  shown  in  Fig. 
3,  by  keeping  the  same  poling  voltage,  we  have  found  that  as  time  progresses,  the  lost  SHG 
intensity  recovers  and  maintains  at  the  same  level  as  before  the  temperature  change.  One 
notes  that  as  the  sample  (5  wt%  NAS)  is  heated  gradually  from  100  to  160°C  (about  50°C 
above  Tg),  SHG  intensity  remains  nearly  constant,  in  contrast  to  the  expected  intensity 
decrease  as  the  temperature  increases.  This  intensity/ temperature  behavior  differs  from  the 
result  of  Mortazavi  et  al,7  who  have  reported  a  strong  temperature  dependence  in  the  main 
chain  polymer.  The  present  result  can  be  understood  as  due  to  electrodes  attraction 
associated  with  the  increase  of  the  effective  poling  field  strength.  As  the  temperature 
increases  above  Tg,  the  compliance  of  the  polymer  increases  and  the  more  compliant 
polymer  facilitates  the  electrodes  separation  to  decrease.  The  decrease  in  the  electrodes 
separation  increases  the  effective  poling  field  strength  which  probably  compensates  the 
O/T)  temperature  dependence  factor  provided  that  a  sufficient  amount  of  time  is  allowed 


Page  7 


for  the  polymer  compliance  to  change, 

2.  Relaxation  of  Polar  Order  Parameter 

Shown  in  Fig.  4A  is  the  decay  curve  of  the  SHG  intensity  for  the  5  wt%  sample  at 
100°C  after  the  poling  field  (=  217  V/^m)  is  removed.  One  notes  the  initial  rapid  intensity 
decrease  at  short  times,  followed  by  a  gradual  decay  at  subsequent  times  in  the  SHG  time 
dependence  curve.  The  nature  of  the  initial  rapid  intensity  decrease  is  not  clearly 
understood,  despite  the  fact  that  it  has  been  attributed  to  the  third  order  nonlinear  optical 
effect8  and  to  the  relaxation  of  the  oriented  chromophores  in  the  liquid-like  region  for  the 
polarized  sample  below  Tg.9  We  have  recently  carried  out  detailed  experiments,  and  have 
found  that  the  portion  of  the  rapid  SHG  intensity  decrease  can  be  made  to  diminish  if  the 
field  higher  than  Es  is  employed  to  pole  the  sample  above  Tg  for  a  sufficiently  long  time. 
In  Fig.  4B,  we  show  the  decay  curve  for  the  sample  poled  at  121°C  using  the  same  poling 
voltage  for  more  than  24  hr.;  one  notes  that  there  is  only  a  negligible  initial  drop  in  the 
SHG  intensity  in  this  case.  The  result  of  the  comprehensive  study  of  the  poling  dynamics  will 
be  published  elsewhere.10 

The  entire  SHG  relaxation  curve,  including  the  initial  intensity  drop,  can  be  fit  to  a 
KWW  (Kahrausch-Williams- Watts)  function:  Ae('(l/T°^f  where  r0  is  the  characteristic 
relaxation  time,  and  p  is  the  parameter  signifying  the  distribution  of  relaxation  times.  By 
fitting  to  one  KWW  function,  one  obtains  the  stretch  exponential  parameter  p  about  the 
order  of  0.3,  which  is  consistent  with  the  previous  result.11  However,  recent  work  has  shown 
that  the  decay  of  the  second  order  nonlinear  susceptibility  is  better  fit  to  a  bi-exponential 
form.12,13  Using  the  birexponential  fit,  we  write  the  SHG  intensity  relaxation  curve  as: 

ISHgW  -  (Ae-,A‘  +  Be'/’2 f  (8) 

where  A  and  B  are  constants  associated  with  the  amplitudes  of  the  initial  and  the  second 
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long  decay,  respectively.  As  mentioned  above,  if  a  sufficiently  high  voltage  is  used  to  pole 
the  sample  above  Tg  for  a  long  time,  the  amplitude  of  the  initial  decay  becomes  rather 
small.  As  a  result,  the  long  decay  associated  with  the  reorientation  of  the  chromophores 
dipoles  is  found  to  be  described  by  a  single  exponential  given  by 

JshgO)  “  r2  exP  (-2t/T 2)  (9) 

The  continuous  curves  shown  in  Fig.  (4A)  and  (4B)  are  the  single  exponential  fit  to 

the  experimental  data.  The  fit  fo  Fig.  (4A)  is  carried  out  by  removing  the  portion 
corresponding  to  the  initial  rapid  intensity  decrease.  One  notes  that  by  excluding  the  initial 
intensity  decrease,  the  single  exponential  function  gives  a  good  fit  to  the  experimental  data. 

The  decay  of  the  SHG  intensity  corresponds  to  the  relaxation  of  the  polar  order  of 
the  NLO  chromophore  dipoles.  The  relaxation  of  the  orientational  order  parameter  can  in 
general  be  expressed  in  terms  of  time  correlation  functions,  but  to  our  knowledge,  this 
connection  has  not  yet  been  shown.  We  now  use  linear  response  theory  to  relate  the 
relaxation  of  the  order  parameter  associated  with  the  decay  of  the  SHG  signal  to  the 
specific  time  correlation  functions. 

To  apply  linear  response  theory,  we  consider  a  collection  of  electric  dipoles  of  NLO 
chromophores  jgA>  y.2  •  *  •  *  #n>  eac^  °f  these  interacts  with  the  external  poling  field  Ep(t) 
and  also  with  the  others.  The  poling  field  depends  on  time  as  it  is  switched  off  after  the 
sample  is  polarized  and  the  relaxation  effect  is  to  be  probed.  The  Hamiltonian  describing 
the  interaction  of  the  dipoles  with  the  poling  field  is  given  by 

Hext(t)  =  -  £  juj  *Ep(t)  =  -  mEp(0  E  cos0j  (10) 


where  0j  is  the  polar  angle  of  dipole  j  with  respect  to  the  poling  field  Ep.  jx  is  the  amplitude 


Page  9 


of  the  electric  dipole  of  the  NLO  chromophore.  In  addition  to  the  external  Hamiltonian,  the 
material  system  is  described  by  a  Hamiltonian  H0  before  the  poling  field  is  applied. 
Following  a  series  of  standard  steps  needed  in  linear  response  theory,14  we  find  that 

<cos”91(0>e=  ♦»(«')  Ep(t')dt'  (11) 

where  the  response  function  0(t,t')  is  given  by  (in  terms  of  classical  mechanics) 

E  <cos”e1(t)  ±  COS  9j  (t')>  (12) 

In  accordance  with  Eqs.  (3)  and  (4),  we  need  03(t,t')  in  order  to  calculate  L^t);  to 
calculate  L^t)  we  need  03(t,f)  and  0j(t,t').  To  derive  Eqs.  (11)  and  (12),  we  have 
assumed  that  initially  the  system  has  no  polar  order.  Polar  order  is  induced  after  the 
external  electric  field  is  turned  on.  However,  as  previously  shown,  under  the  weak  poling 
field  conditions  only  the  <P1(cos0)>  type  of  polar  order  is  excited.  The  excitation  of  the 
<P3(cos0>  order  requires  a  high  poling  field  and  cannot  be  calculated  with  linear  response 
theory.  We  shall  in  another  publication  consider  the  high  field  case  and  consider  here  only 
the  low  poling  field  condition. 

We  now  substitute  Eq.  (12)  into  Eq.  (11)  (with  n  =  1)  and  carry  out  the  integration 
by  part  to  reduce  Eq.  (11)  to 
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(13) 


N 

<  cosS^t)  >E  =  JL  £  <  cos0j(t)  cos0j(t)  >  Ep(t) 
kT  jml 

■if  f  <cos91(t)  cos9:(t-)>  ep(f)df 
kX  *  *• 

where  the  dot  above  Ep(t')  designates  taking  the  derivative  with  respect  to  t'. 

We  next  consider  the  situation  that  after  the  NLO  polymer  material  is  poled  for  a 
sufficiently  long  time,  the  dipolar  order  has  reached  a  steady  state.  At  t=0,  we  turn  off  the 
poling  field  and  begin  to  investigate  the  relaxation  behavior  of  the  established  polar 
orientational  order.  The  investigation  is  proceeded  by  observing  the  decay  of  the  SHG  signal 
at  t>0.  Namely,  we  set 

EP(t)  =  Ep  0(t)  (14) 

where  0(t)  equals  1  for  t<0,  and  equals  0  for  t>0.  From  Eq.  (12),  we  obtain 

Ep(t)  =  -  Ep<S(t)  (15) 

where  <S(t)  is  the  Dirac  delta  function.  Substituting  Eqs.  (14)  and  (15)  into  Eq.  (13),  we 
obtain  for  t>0 


<cos01(t)>E  =  (-^)  <cos®i(t)  cos  0j(O)> 

=  a  [  < cos0^(t)  cos01>  +  (N-l)  <  cos0j(t)  cos  02>] 


(16) 


where  the  angular  brackets  without  the  subscript  E  indicate  that  the  ensemble  average  is 
taken  without  the  poling  field.  The  first  term  on  the  right  hand  side  of  Eq.  (13)  is  absent 
for  t>0  as  the  poling  field  is  removed  at  t=0.  In  Eq.  (16),  we  have  succeeded  in  relating  the 
decay  of  the  polar  order  parameter  to  time  correlation  functions.  One  notes  in  Eq.  (16) 
there  are  two  types  of  time  correlation  functions,  one  corresponding  to  the  self  time 
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correlation  and  the  other  pair  time  correlation.  In  the  previous  work,  only  the  analysis  of 
the  self  time  correlation  function  is  considered.15  We  emphasize  in  this  paper  the 
importance  of  the  pair  correlation  contribution  to  the  decay  of  the  SHG  intensity.  Without 
the  pair  correlation  term,  it  would  be  impossible  to  understand  the  concentration 
dependence16  and  the  present  poling  field  dependence  of  the  relaxation  time. 


3.  Effect  of  Pair  Correlation 

In  Eq.  (16),  we  have  established  that  the  time  dependence  of  the  order  parameter 
<  cosS^t)  >E  is  related  to  the  correlation  functions  <cos0j(t)  cos01>  and 
<cos0l(t)  cos02>.  To  show  the  importance  of  pair  correlation  contribution,  we  consider  a 
model  with  two  interacting  dipoles  immersed  in  the  environment  of  the  amorphous  polymer 
above  its  glass  transition  temperature. 

The  dipoles  are  initially  oriented  by  the  poling  field.  After  the  poling  field  is 
removed,  the  dipoles  undergo  reorientational  motion  to  randomize  the  orientational  order 
established  by  the  poling  field.  The  reorientational  dynamics  of  the  dipoles  can  be  described 
by  the  generalized  Langevin  equation.17 

If  we  denote  the  state  of  dipoles  by  the  dynamic  variables  Ql  and  Q2,  then  the 
kinetic  equation  for  the  vector  variable 


A  = 


(17) 


is  given  by  the  generalized  Langevin  equation17  (in  the  Markov  approximation) 
where  the  relaxation  matrix  is  given  by 
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|A(t)  =  r-A(t)+F(t) 

&  (18) 

(Kll  K12 

“W  *22, 

Here  Ky  are  the  Onsager  kinetic  coefficients  associated  with  the  random  force  F(t).  Since 
the  two  dipoles  are  assumed  to  be  chemically  equivalent,  we  set  Kn  =  K22  =  d,  and  K12 
=  K21  =  h.  Here  d  is  associated  with  the  decay  rate  of  uncorrelated  single  dipole,  and  h 
represents  dynamic  correlation  of  the  two  dipoles.  In  addition,  Qx  and  Q2  are  proportional 
to  cosOj^  and  cos02,  respectively,  we  set  <  |Q1|2>  =  <  |Q2|2>  and 
<QiQ2*>  =  <Q2Q-|*>  =p.  Here  p  represents  the  static  pair  correlation  between  the  two 
dipoles.  Using  Eq.  (6),  we  can  write  for  the  case  of  two  dipoles  the  initial  orientational 
order  parameter  before  turning  off  the  poling  field  is 

a  +  p  =  LjjCO)  =  [<cos201>  +  <cos01cos02>]  (20) 


<|Q|I>  <QlQ2)l 
<Q,Q>  <!Q,I2> 


(19) 


Equation  (20)  will  be  used  as  the  initial  condition  for  the  kinetic  equation  given  by 
Eq.  (18).  Thus,  in  order  to  show  that  Eq.  (18)  can  be  used  to  describe  the  decay  of  the  polar 
orientational  order  parameter  after  the  poling  field  is  turned  off,  we  need  to  make 
connection  of  Eq.  (18)  to  Eq.  (16). 

We  take  the  time  derivative  on  both  sides  of  Eq.  (16)  and  obtain  after  associating 
Qj  with  cos0,  (i  =  1  and  2): 
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~  <COS01(t)>E 
ot 


a  d 

- -  (<COS01(t)  COS01>  +  <COS0j(t)  COS02> 

2  dt 


+  <COS02(t)  COS01>  +  <COS02(t)  COS02>] 

=  14  [<Q1(t)Q1>  +  <Q2(t)Q2>  +  <Q2(i)Qi>  +  <Q2(t)Q2>]  (21) 

2  dt 


where  we  have  added  two  terms  associated  with  the  representative  dipole  2  to  exploit  the 
symmetry  of  the  problem.  The  initial  condition  for  Eq  (21)  is  given  by  Eq.  (20). 

The  calculation  of  the  sum  of  four  correlation  functions  on  the  right  hand  side  of  Eq 
(21)  is  best  proceeded  by  considering  the  correlation  matrix 


C(.)  = 


'<Q,(t)Q,>  <Qi(t)Q3>' 
<Qj(t)Q,>  <Qj(0Qi> 


(22) 


The  equation  of  motion  for  £(t)  is  obtained  from  Eq  (18)  by  multiplying  the  adjoin  of  A 
to  the  right  side  of  both  equations  and  then  taking  the  ensemble  average: 


4  ff®  ■  -£  -gw 


(23) 


Unlike  Eq.  (18),  the  random  force  F  (t)  is  absent  explicitly  in  Eq.  (23)  owing  to  the  fact  that 
the  correlation  between  random  force  F  and  dynamic  variable  A  vanishes.17 

The  sura  of  the  diagonal  elements  of  £(t),  [<Qj(t)  Qt>  +  <Q2(t)  Q2>],  and  the 
sum  of  off-diagonal  terms  of  £(t),  [  <  Qt(t)  Q2>  +  <  Q2(t)  Qx  >  j,  are  known  as  the  self-time 
correlation  functions  and  the  pair-time  correlation  functions,  respectively.  The  orientational 
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order  parameter  can  thus  be  obtained  from  the  elements  of  the  time  correlation  function 
matrix. 


To  calculate  the  time  correlation  function  matrix,  we  take  the  Laplace  transform  of 
Eq  (23)  and  convert  it  into  an  algebraic  equation: 


(24) 


where  X  is  a  2  by  2  unit  matrix,  and 


Coo  -  /■  dte'“  C(t) 


(25) 


£(0)  is  £(t)  at  t  =  0.  Using  the  definition  of  thej  matrix  (Eq  (19)),  we  have  found  after 
performing  several  straightforward  matrix  multiplication  steps  that 


Coo  - 


1 


(z+L)2-M2 


a(z+L)-pM  p(z+L)-aM 
[p(z+L)-aM  a(z+L)-pM 


(26) 


where 


L  = 


a2p2 


(ad  -  ph) 


(27  A) 


M  = 


a2p2 


(ah  -  pd) 


The  self-  and  pair  correlation  parts  of  Eq  (26)  are  thus  given  by 


(27B) 


C.O0  -  {a  (z  +L)  -pM} 

(z  +L)2  -M2 


(28A) 
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(28B) 


One  notes  from  Eq  (16)  that  the  Laplace  transform  of  the  order  parameter  <  cosS^t)  >E 
is  equal  to  one  half  of  the  sum  of  Cs(z)  and  Cp(z);  i.e., 


<cos0t(z)>B = dte  '*  <cos91(t)>B  =  ~  [6,(z)  +  Cp  (z)] 

=  «  +  P  =  («  +  P) 

z+L+M  z  +  (d+h)/(a +  P) 


(29) 


Equation  (29)  can  be  reverted  to  time  by  taking  the  inverse  Laplace  transform.  Thus,  we 
have  found  that  the  decay  of  the  POP  to  the  polar  orientational  order  parameter  is  given 
by 

<  cosQ^t)  >E  =  (a  +  j8)  e'*/T  (30) 

where  the  relaxation  r  is  given  by 


T  -  («  +  P> 

(d+h) 


(31) 


Equation  (31)  shows  that  the  relaxation  time  r  for  a  two-coupled  dipole  system  is 
proportional  to  a  +  p,  which  in  accordance  with  Eq  (20),  is  equal  to  the  order  parameter 
1^(0).  Thus,  the  relaxation  time  increases  as  the  polar  orientational  order  parameter 
increases.  Both  kinetic  coefficients  d  and  h  also  contribute  to  the  relaxation  time  r,  but  they 
do  not  depend  on  the  order  parameter.  It  is  important  to  point  out  that  if  the  pair 
correlation  contribution  (Eq.  (28B))  is  not  included  in  the  calculation,  the  decay  of  the  POP 
will  only  be  given  by  the  self-part  (Eq.  (28A)),  which  does  not  have  the  simple  form  of  Eq. 
(30). 


Although  the  above  calculation  is  made  on  a  two-coupled  dipole  model,  we  have 
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shown  elsewhere,18  that  the  above  result  is  only  slightly  modified  when  the  calculation  is 
extended  to  a  general  N-coupled  dipolar  system.  The  modification  corresponds  to  the 
replacement  of  (3  by  (N-l)/3,  and  h  by  (N-l)h. 

To  verify  the  prediction  of  Hq.  (33),  we  have  carried  out  the  measurement  of  the 
SHG  signal  as  a  function  of  time  using  a  sample  (containing  5  wt%  NAS)  under  an 
isothermal  condition.  The  sample  isothermally  controlled  at  U5°C  are  repeatedly  poled  by 
an  electric  field  of  variable  strengths.  The  shape  of  the  decay  curve  at  each  poling  field  is 
fit  to  Eq.  (9)  to  obtain  the  relaxation  time  r.  In  Fig.  5,  we  plot  the  relaxation  time  r  versus 
the  strength  of  the  poling  field  Ep.  One  notes  that  r  increases  with  Ep,  consistent  with  the 
prediction  of  Eq.  (31).  As  a  +  j3  is  equal  to  the  order  parameter  L31(0),  which  increases 
with  the  poling  field  strength. 

We  have  thus  shown  that  the  decay  of  the  SHG  signal  can  be  impeded  by  increasing 
the  poling  field  strength  Ep.  The  increase  in  Ep  results  in  a  larger  POP,  which  thus  stabilizes 
the  SHG  signal. 

In  summary,  we  have  investigated  the  electric  field  dependence  of  the  SHG  intensity 
and  its  stability  using  a  guest/host  polymer  system.  We  have  shown  that  under  the  low 
poling  field  condition  the  SHG  intensity  is  proportional  to  Ep2,  provided  that  the  poling  field 
is  greater  than  the  threshold  field  Es.  Below  the  threshold  field,  only  the  dipoles  near  the 
electrodes  are  oriented.  The  macroscopic  orientational  polar  order  develops  as  the  poling 
field  exceeds  Es.  The  decay  of  the  SHG  signal  is  strongly  affected  by  the  pair  orientational 
correlation  which  significantly  contributes  to  the  polar  order.  The  lengthening  of  the 
orientational  relaxation  time  by  the  increase  of  the  poling  field  strength  is  observed. 
Theoretical  analysis  utilizing  linear  response  theory  shows  that  the  relaxation  time  associated 
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with  the  decay  of  polar  order  is  proportional  to  the  polar  order  parameter  initially  induced 
in  the  sample.  We  have  carried  out  experimental  measurements  and  verified  the  prediction. 
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Figure  Captions 


Fig.  1  Variation  of  the  glass  transition  temperature  (Tg)  of  the  NAS/PMMA  system  with 
the  NAS  concentration. 

Fig.  2  Amplitude  of  the  SHG  signal  for  the  5  wt%  NAS  in  the  PMMA  sample  at  110°C 
plotted  as  a  function  of  the  square  of  the  poling  field  strength  (Ep2). 

Fig.  3  Temperature  dependence  of  the  SHG  intensity  for  the  5  wt%  sample. 

Fig.  4  The  SHG  decay  data  for  the  sample  poled  at  two  temperatures,  at  the  same  electric 
field.  Note  the  disappearance  of  the  initial  intensity  decrease  in  Curve  B.  The  solid 
curves  are  the  single  exponential  fit  to  the  experimental  data. 

Fig.  5  Relaxation  time  constant  r  associated  with  the  decay  of  the  polar  order  for  the  5 
wt%  NAS  sample  isothermally  controlled  at  115°C  plotted  as  a  function  of  the 
strength  of  the  poling  field.  The  solid  curve  is  drawn  to  indicate  the  trend. 
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